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We find that a system of particles interacting through a simple isotropic potential with a softened 
core is able to exhibit a rich phase behavior including: a liquid-liquid phase transition in the super- 
cooled phase, as has been suggested for water; a gas- liquid- liquid triple point; a freezing line with 
anomalous reentrant behavior. The essential ingredient leading to these features resides in that the 
potential investigated gives origin to two effective core radii. 

PACS numbers: 61.20.Gy, 64.70.-p 
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Fluid-fluid or liquid-liquid transitions have been pro- 
posed for many different substances such as supercooled 
water, C,S,Ga,Se,Te,/2,Cs, S1O2 etc. The transition 
may occur between stable phases as in carbon |||| or 
between metastable phases as has been suggested in su- 
percooled water In general, these substances also 
exhibit different amorphous states as well as multiple 
crystalline structures [Q. The complex intermolecular 
interactions characterizing the above materials are usu- 
ally modelled through anisotropic potentials depending 
on the molecular orientation and thus the possibility 
that simple fluids interacting through isotropic potentials 
may exhibit similar exotic phase behaviors is a challeng- 
ing field of investigation. 

As argued by Mishima and Stanley |I| the competi- 
tion between expanded and compact structures in fluid- 
fluid (liquid-liquid) transitions suggests that the poten- 
tial should possess two equilibrium positions. The most 
obvious form with such feature is one with two wells. 
Such potentials were shown to give rise to water-like ther- 
modynamic anomalies, though the presence of a new crit- 
ical point could not be directly observed ||. Another 
form of interparticlc interaction which could produce dif- 
ferent equilibrium positions is that in which there is a 
region of negative curvature in the repulsive core: these 
so-called softened-core potentials were proposed by Stell 
and Hemmer Q who argued that they might produce 
a second transition if a first already exists. Recently, 
through a mixed numerical-mean field type calculation, 
it was found that a potential consisting of a softened-core 
plus an infinite range van der Waals attractive term may 
give rise to a second critical point Very recently, 

molecular dynamics simulation showed for a softened- 
core potential with an attractive square well evidences of 
a transition between two fluid phases in the supercooled 
region p2| . 

The purpose of this Letter is to report the findings of a 
study of the phase behavior of a system of particles inter- 
acting through a potential consisting of a hard core with 
a finite repulsive shoulder and an attractive square well. 
Our analysis, based mainly on thermodynamically self- 
consistent (TSC) integral equations for fluids pf and 



partly on Monte-Carlo (MC) simulations, shows the ex- 
istence of a liquid-gas critical point in the stable fluid 
phase and of a liquid-liquid critical point in the super- 
cooled region. The liquid-gas and liquid-liquid coexis- 
tence lines meet in a gas- liquid-liquid triple point. More- 
over, the behavior of the freezing line, estimated through 
one-phase criteria [jl3|, such as the Hansen- Verlet (HV) 
rule |l6| and the entropic criterion based on the anal- 
ysis of residual multiparticle entropy p7| , is consistent 
with the existence of multiple crystalline structures in 
the solid phase. This is the first time that a microscopic 
theory linking directly JIJ] the behavior of the system to 
the form of the interparticle pair potential, predicts for a 
simple fluid the existence of a liquid-liquid critical point 
and of a gas-liquid- liquid triple point. 

The chosen potential has a repulsive part V rep (r) con- 
sisting of a hard core of radius ro — a and a repulsive 
square shoulder of height e and radius r\ — 2.5a, plus 
an attractive component V a tt r {r) having the form of a 
square well of depth 1.25e extending from r\ = 2.5<r to 
r2 = 3ct. To reach a thorough comprehension of the role 
played by the different components of the potential we 
first study its purely repulsive part, and then consider 
the effect of adding the attractive component. 

As a preliminary investigation we study through MC 
simulation M ] a system of particles interacting through 
the potential V rep (r) in 2D, since in this case the spa- 
tial arrangement of the particles can be easily visualized. 
Snapshots representing the configuration of the system 
at a low temperature and different densities are shown in 
Fig.l, together with the corresponding structure factors 
S(k). As the density increases, the system rapidly under- 
goes a very strong ordering (Fig. 1(b)) indicated by the 
huge first peak of S(k) (the splitted second peak of S(k) 
can be related to the presence of lattice dislocations evi- 
dent in the figure). Upon further increasing the density, 
particles begin to enter the soft core and many hard cores 
get in contact: the first peak of S(k) decreases while the 
second peak of S(k) is considerably enhanced (Fig. 1(d)). 
Finally at very high densities (not shown in the figure) 
the third peak of S(k) becomes the highest one, while 
the first two peaks are greatly depressed. 
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The progressive "rising and falling" of the peaks of the 
structure factors reflects the turning on/off of different 
effective length scales. When the temperature T and the 
density p are sufficiently small, the soft core is practically 
impenetrable and the particles behave as hard spheres of 
radius r%. As T and p increase, more and more parti- 
cles penetrate the soft core until this becomes scarcely 
influent and the system is essentially equivalent to an 
assembly of hard spheres of radius r®. In general, the 
system can be considered a "mixture" of two popula- 
tions of hard spheres, of radius r$ and r\ respectively. 
The relative concentration of the two species is fixed by 
the values of T and p. Thus, in contrast to standard 
simple fluids, the system has three possible length scales: 
r ii r io = ( r i + r o)/2 and ro, and as many indicators of 
structural ordering, namely the peaks of the structure 
factor corresponding to the wavevectors k\, k\$, ho, as- 
sociated to these lengths. We stress the crucial role of 
the finite repulsive shoulder; due to its presence particles 
can be in one of two "states" : this is the essential feature 
which opens the possibility of liquid-liquid-immiscibility 
in a pure substance. 

Let us now consider a system of particles interacting 
through the potential V rep (r) in 31?. We study its struc- 
tural and thermodynamical properties using the TSC 
Roger- Young (RY) integral equation p3[ . According to 
HV rule a fluid is expected to undergo crystallization 
when the first (main) peak of S(k) attains the value 2.85 
p6[ . This statement refers to simple fluids with a single 
length scale. In our case, different length scales come into 
play and one must consider all the associated indicators 
of structural ordering. As shown in Fig. 2, their behav- 
ior, upon increasing the density at constant T, appears 
analogous to that observed in 2D, with each peak in turn 
rising and then "falling" down (except that associated to 
the hard core radius which becomes higher and higher 
as p increases). The overall trend is confirmed by MC 
calculation, though the theory j24| tends to overestimate 
the heights of S(k\o) and S(k ) at high densities. 

In Fig. 3 we show the loeii of the points of the plane T,p 
for which S(ki), 5(fcio) and S(ko) are equal to 2.85. The 
freezing line predicted through a straightforward appli- 
cation of HV rule to the "anomalous" simple fluid inves- 
tigated, should coincide with the line which bounds the 
region where at least one of the peaks of the structure 
factor is greater than 2.85. This line shows a reentrant 
behavior which suggests that the system may undergo 
crystallization upon the density decreasing. The behav- 
ior can be related to the presence of the finite repulsive 
shoulder. In fact, the fraction of particles penetrating the 
soft core increases with the density (at constant T) thus 
"generating" additional space for the system. This effect 
is particularly important where one length scale begins 
to become less effective in favor of the smaller one. In 
these regions the phenomenon may overcompensate the 
general decrease of the space available to the system upon 



the density increasing thus causing a tendency of the sys- 
tem to become less ordered: accordingly, the freezing line 
may have a negative derivative. We note, however, that 
this behavior might be overestimated in Fig. 3. In fact, 
since different length scales are contemporarily effective 
and contribute to the global structural ordering of the 
system, we expect that a more realistic freezing line has 
a smoother behavior which "blends" together the upper 
portions of the different curves of Fig. 3. 

The inset in Fig. 3 shows a magnification of the low 
temperature-low density region. Here, since only one 
length scale (t\) is effective, the prediction of HV is ex- 
pected to be reliable. The freezing line starts nearly ver- 
tical at a density p* ~ 0.06 which corresponds to the 
freezing density pr\ ~ 0.945 of a fluid of hard spheres of 
radius t\. As p increases, the freezing line bends and 
shows a reentrant behavior: the phenomenon is asso- 
ciated with the very onset of the soft core penetration 
which, for the reasons above discussed, has a disordering 
effect on the system. This is reflected in the anoma- 
lous behavior shown by the isothermal compressibility 
Xt, which suddenly increases with the density (see same 
figure). Correspondingly, the first peak of S(k) abruptly 
decreases: this occurs also at lower temperatures (where 
the value of S(ki) remains however greater than 2.85) 
along a line coinciding approximately with the extension 
of the portion of the freezing line with negative derivative 
p6[ . This line meets the T = axis at a density p* ~ 0.09 
which corresponds to the closest packing of hard spheres 
of radius r\ (occurring at pr\ — \/2). The above results 
lead to conclude that the region shadowed in the inset 
of Fig. 3 corresponds to an expanded solid phase of the 
system. 

We now investigate the phase behavior of a system 
of particles interacting through the full potential V rep (r) 
plus V a ttr{r)- Calculations a re p erformed making use of 
the HMSA TSC equation |j, which is better suited 
than RY equation for interparticle interactions including 
an attractive component and reduces to this for purely 
repulsive potentials. The phase diagram of the system 
is shown in Fig. 4. Two coexistence curves occur [p8| , 
each terminating at a critical point, denoted CI and C2. 
The critical densities and temperatures are respectively 
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0.06, T. 



ci 



1.3, and p* C2 = 0.77, T* 2 = 0.55. 



These values were estimated using the rectilinear diam- 
eter rule and the scaling relationship for the width of 
the coexistence curve with the non-classical exponent 
j3 « 0.325 |29]. Below Tq\ the system separates into 
a gas and a liquid phase. The liquid phase is not unique 
since, below Tci {Tc2 < Tci), separates into distinct 
low-density (LD) and high-density (HD) phases. Since 
the critical point C2 is well below the freezing line, the 
liquid-liquid transition occurs between metastable phases 
in the supercooled region of the system. This feature re- 
calls the scenario proposed for water E 19 1, but in that 
case the liquid-liquid coexistence line is expected to start 
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from C2 running at higher pressures as T decreases ||. 
In the system investigated, the contrary is observed, this 
line running at lower pressures as T decreases (see inset 
of Fig. 4). This makes possible a new feature: the simul- 
taneous coexistence of three fluid phases. In fact, the 
gas-liquid and the liquid-liquid coexistence lines meet in 
a gas-liquid-liquid triple (GLL) point which lays in the 
supercooled phase (pg LL ~ 0.57, Tq LL 0.53). 

The addition of the attractive well causes a shift to- 
wards higher temperatures of the locii of the points where 
S(ki), S(kio) and S(ko) are equal to 2.85, while their lo- 
cation in density remains essentially unaltered. Only a 
small portion of the line S(k\) = 2.85 is visible since it 
lays almost entirely in a region, correspondingly approx- 
imately to the gas-liquid spinodal decomposition, where 
the theory is unstable. The freezing line corresponding to 
HV rule and that estimated through the entropic crite- 
rion are in substantial agreement at low and high densi- 
ties (pressures) , yielding close estimates of the gas- liquid- 
solid {GLS) triple point {p* GLS w 0.29, T^ LS w 0.97). 
On the contrary, these lines differ markedely in the inter- 
mediate region, where the first one shows a very evident 
reentrant behavior, while in the second approach this can 
be appreciated only numerically . For the reasons ex- 
posed in the purely repulsive case, the second result is 
probably more reliable. Thus, as shown in the inset in 
Fig. 4, the freezing temperature increases initially with 
pressure, then remains essentially constant (according to 
the entropic criterion) in the pressure range 1 < P* < 2, 
(P* = Pa 3 / e) and eventually increases again with pres- 
sure. Though our estimate of the freezing line is based 
solely on one-phase criteria, its shape, with branches hav- 
ing distinctly different slopes, is consistent with the pos- 
sibility that structural changes occur in the solid state 
of the system. Consequently, transitions may be possible 
between solid phases of the system investigated. 

The results presented in this Letter show that a pure 
model system, with a softened-core isotropic potential, 
may have a rich phase behavior with features typi- 
cal of substances characterized by much more complex 
anisotropic interactions. The possibility that simple sub- 
stances related to the model investigated may exist (or 
may be "realized" ) in nature well deserves to be investi- 
gated. 
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aquinta, C.Caccamo, S.Dugdale, for helpful suggestions 
and for interesting and stimulating discussions. One of 
us (G.M.) wishes to thank H.E.Stanley for stimulating 
his interest on the subject. 
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FIG. 2. Purely repulsive potential. Structure factors within the RY equation and MC simulation. 
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FIG. 3. Purely repulsive potential. The dashed lines are the locii of points of the T,p plane where S(ki) (full triangles), 
S(feio) (squares) and S(ko) (open triangles) are equal to 2.85. The inset shows a magnification of the low T-low p region; the 
isothermal compressibility \t (thick solid line) is calculated along the isotherm T* = 0.103 (xo is the ideal gas compressibility). 
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FIG. 4. Phase diagram in the T,p plane. Coexistence lines: gas- liquid (circles), liquid-liquid (dots). The full diamonds 
represent the critical points. The dashed lines are the locii of points where S(ki) (full triangles), S(fcio) (squares) and S(ko) 
(open triangles) are equal to 2.85. The solid line with no symbols is the freezing line estimated through the entropic criterion. 
Inset: P,T phase diagram; pressure is given in units of e/a 3 . 
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